home *** CD-ROM | disk | FTP | other *** search
/ Turnbull China Bikeride / Turnbull China Bikeride - Disc 1.iso / ARGONET / PD / MATHS / RLAB / RLAB125.ZIP / !RLaB / examples / 3d next >
Text File  |  1995-05-20  |  1KB  |  69 lines
  1. plstart();
  2.  
  3. # Make 3D data
  4.  
  5. # Sombrero
  6.  
  7. NX = 20;
  8. NY = 20;
  9. xx = zeros (1, NX);
  10. yy = zeros (1, NY);
  11. for (i in 1:NX) { xx[i] = (i - NX/2)/(NX/2); }
  12. for (i in 1:NY) { yy[i] = (i - NY/2)/(NY/2); }
  13.  
  14. zz = zeros (NX, NY);
  15.  
  16. for (i in 1:NX)
  17. {
  18.   for (j in 1:NY)
  19.   {
  20.     r = sqrt (xx[i]^2 + yy[j]^2);
  21.     zz[i;j] = exp (-r * r) * cos (2*pi*r);
  22.   }
  23. }
  24.  
  25. # Sin - Cos surface
  26.  
  27. x1 = -3:3:.2;
  28. y1 = -3:3:.2;
  29. z1 = zeros (x1.n, y1.n);
  30.  
  31. for (i in 1:x1.n)
  32. {
  33.   for(j in 1:y1.n)
  34.   {
  35.     z1[i;j] = sin(y1[j]) * cos(x1[i]);
  36.   }
  37. }
  38.  
  39. # Slanted plane
  40.  
  41. x2 = 1:10;
  42. y2 = 1:10;
  43. z2 = zeros (x2.n, y2.n);
  44. for (i in 1:x2.n) { z2[i;] = i*ones(1,y2.n); }
  45.  
  46. # Now create some plots
  47.  
  48. pltitle("PLOT3");
  49. plot3 (<< x = xx; y = yy; z = zz>>);
  50. pause();
  51. pltitle("PLMESH");
  52. plmesh (<< x = xx; y = yy; z = zz>>);
  53. pause();
  54.  
  55. pltitle("PLOT3");
  56. plot3 (<< x = x1; y = y1; z = z1>>);
  57. pause();
  58. pltitle("PLMESH");
  59. plmesh (<< x = x1; y = y1; z = z1>>);
  60. pause();
  61.  
  62. pltitle("PLOT3");
  63. plot3 (<< x = x2; y = y2; z = z2>>);
  64. pause();
  65. pltitle("PLMESH");
  66. plmesh (<< x = x2; y = y2; z = z2>>);
  67. pause();
  68.  
  69.